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Abstract 

We propose a mechanism to have a smooth transition from a pre-Big Bang phase to 
a standard cosmological phase. Such transition is driven by gravitational production of 
heavy massive string states that backreact on the geometry to stop the growth of the 
curvature. Close to the string scale, particle creation can become effective because the 
string phase space compensate the exponential suppression of the particle production. 
Numerical solutions for the evolution of the Universe with this source are presented. 



1 Introduction 



The standard cosmological model, even in its inflationary extension, is plagued by an "ini- 
tial condition problem" [1], since it must emerge from an initial singularity. One of the 
most successful proposal to overcome this problem was the introduction of the so-called 
pre-Big Bang scenario [2] (for a review, see [3]), developed in the context of string theory, 
in which the Universe is supposed to emerge from a string vacuum (so the initial geom- 
etry is flat), perturbed by quantum fluctuations of the dilaton (other proposals were put 
forward in different contexts, such as brane cosmology and Loop Quantum Gravity, see 
for example [4-10]). This is achieved by implementing purely stringy symmetries on the 
string cosmology equations. Those symmetries suggest the existence of a phase of growing 
curvature. However, it is not clear, a priori how to connect smoothly the pre-big bang and 
the post-big bang phases, since, at the classical level, the two phases are still disconnected 
by a singularity. Many ideas were proposed to overcome this problem, all involving further 
assumptions on the dynamics of the Universe, namely higher order correction (either loop or 
a' contribution) and negative energy density contribution to obtain the desired "repulsive 
gravity" (an incomplete list would include [11-16]). 

In this paper we will propose a mechanism to obtain a regular evolution of the curvature 
that rely on production of heavy massive string states by gravitational backreaction. It is 
well known [17] that an evolving Universe will produce particles because of the squeezing 
of the ground state. It has been speculated in the literature that this effect, though being 
negligible, could actually have some physical consequences [18-20]. The idea is that the 
exponential suppression of the energy density of particles created during the evolution of 
the Universe is compensated by the exponential growing of the multiplicity of states above 
the string scale (the Hagedorn spectrum), thus giving a sensible contribution to the total 
energy density. Of course, dealing with gravitationally induced particle production brings 
one to deal with a sort of semiclassical realization of quantum gravity, which is not under 
control, and in fact the authors themselves of the mentioned papers admit some of the 
ideas are rather speculative. Here we want to go one step forward with speculations, and 
assume that the gravitationally produced energy backreacts locally on the geometry of the 
Universe. We will see that even a rough estimate such as the one we are going to show here 
is enough to stop the growing of the Ricci curvature and induce a bell-shaped behaviour. 
Let us stress that all the analysis is performed under the assumtpion that the universe stays 
in a low-energy (low curvature) regime in which supergravity is a valid approximation, and 
the only ingredient taken from stringy physics is the behaviour of the high energy spectrum 
that influences the phase space of low curvature universe. In fact, we are going to treat 
string massive modes as scalars, assuming however that their multiplicity is controlled by 
the Hagedorn spectrum. The validity of this assumption is verified by checking that the 
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curvature found vith the modified cosmological equations do not exceed the string scale 
(which has been set to unity in the numerical evaluations!. 

We are aware that some of the passages we will show, though physically well-motivated, 
can be challenged on a formal basis. Nevertheless, since our goal was only to check if 
a backreaction mechanism could provide for a graceful exit without invoking any further 
assumptions and higher order modifications of the string cosmology equations, we just 
skipped over formal complexities to have an estimate of the order of magnitude of the 
effect. We plan to come back on all the mathematical issues (among other things we will 
discuss below) in a forthcoming paper. 

The paper is organized as follows: Section 2 is devoted to some formal developments that 
will lead us to write the enregy density (and pressure) produced by gravitationally excited 
string states in term of a time dependent Bogoliubov coefficient. In section ?? we will use 
this expression to find, under certain assumption, a sound expression for the energy density, 
to be put into the string cosmology equations whcih are solved numerically in section 4. 
Finally, in section 5 we will comment on our results and the (rather strong) approximations 
assumed in obtaining them, and on how to further develop the present study in several 
directions. 

2 QFT approach to gravitational backreaction 

Let us then consider a massive scalar field living in a curved space. The action is 

S = \ J d 4 x^/\g~\ (g^d^dA - mV) • (1) 

Let us assume that the geometry of the Universe is a conformally flat FRW; so, in the 
conformal time gauge, the metric is 

ds 2 =a 2 (r]) (dri 2 - dx 2 ) , (2) 

and the action can be written as 

S = I J d V d 3 x a 2 {rj) (cf)' 2 - (Vcf>) 2 - m 2 a 2 (r))(f) 2 ) , (3) 

where (V(/>) 2 stands for Ylii^i^) 2 - Now we introduce the canonical field ip(rj, x) = a(ry)0(ry, x), 
so that the action can be rewritten as 



S = ^ / dr]d 3 x 



(4) 



This action can be interpreted as describing a scalar field living on a flat space-time, with a 
potential which controls its gravitational interactions. From the action (4) we can evaluate 
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the (canonical) energy momentum tensor. 



a" 



{d^f - ( mV - - ) ^ 2 



a 



V' (5) 



The diagonal elements of this tensor are the energy density and the pressure. The field 
conjugate to ip is 

so the Hamiltonian density, to be identified with the energy density, and the pressure along 
the i th direction are, according to (5) 

Pm(r/,x) = - [vr 2 + (Vc^) 2 + V^cp 2 ] , 
Pm,i(ri,x) = \ [** + (y^f ~ V(r,)ip 2 ] , (7) 

where, of course, V(rj) = m 2 a 2 (rj) — 

Now we expand the fields in Fourier modes and, following conventional quantization 
procedure, we promote (r/-dependent) Fourier coefficients to operators: 

<^ x ) = J^fj2 J d3k Mv)^, 

TTfax) = j^j~ 2 J d 3 k7T k ( V y^, (8) 

where nk(rj) = (p' k (r]) and ipk(f]) satisfies the (operatorial) evolution equation: 

^(r])+^(v)Vk(r]) = 0, (9) 

with 0Jk(v) = \A 2 + V(rj). Moreover, since the field <p is real, and thus Hermitean, the 
modes must satisfy the relation 

<Pk(v) = <P-k(r)), 

= ( 10 ) 

We can then write down the energy density and pressure in terms of these modes. By 
defining the vector 

&w = f ) , (id 

the densities (7) can be rewritten as 

PmM*) = 5 J -fffirdW Q ^ Cp(r/K (k - p) ' x , (12) 
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with 



k-p + V(ri) 
1 



- ( k ' p - v{ ' ,) [)■ ('»> 

This expression is useful, because we know [21,22] that the evolution for (k(v) is given by: 

(k{v) = U k {v,Vo)(k(vo), (14) 
where the propagator can be written as 

U k ( V ,r,o) = ( Mri,m) ^'^'V (15) 
V C k (v,Vo) Dk{v,Vo) J 

and the four coefficients of the matrix are 

A k (v,Vo) = i [9k(Vo)fk(v) ~ 9* k {vo)fk{v)} , 

B k {v,Vo) = i[fk{v)fk(Vo) ~ fk(v)fk(vo)], 

C k (r), r) ) = i [9k{vo)gUv) ~ 9k(Vo)9k(ri)} , 

D k {v,Vo) = i[9k{v)fk(Vo) ~ 9*kiv)fk{vo)}- (16) 

Here fk(v) 1S solution of the canonical evolution equation (9) for modes of the field tp and 
gk(v) is solution of the mode equation (which we have not reported) for the conjugate 
momentum ir (so, of course, the relation gk(rj) = fki 7 ]) holds), while rjo is a suitable time on 
which we will define initial condition. We will assume that rjo is far enough in the past so 
that it is possible to consider the space to be Minkowski. Using eq. (14) we can thus write: 

p^x) = \J f^d( V0 ) SkMm) C P (%K (k - p) ' x , 

JWfox) = \J ^|rClfo>) rW(i7,i») C P (vo)e^-^, (17) 

where, 

Sk, P {v,Vo) = ul(7],7] )Pk iP Uk{v,Vo), 

T kl(v,vo) = ul(^)Qfjj k (^m)- (18) 

Now we want to express this in terms of creation-annihilation operators. Notice that on an 
actual Minkowski space, there are no ambiguities in building a Fock space; so we define: 

/ l \ / l \ 



^>=(4^hHvtr L *' (19) 
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with a k and a\ satisfying the usual canonical commutation rules and u> k = w&(?7o)- Inserting 
this expression in (17) we get 



Pm 
Pm,i 



d 3 k d 3 p 

(2vr) 3 
d 3 k d 3 ° 



Wo + <?( 2 )„t„t , o(3) , o(4) t 

b k,p a k a P + b k,p a k a -p + b k,p a -k a P + b k,p a -k a ~p 



i(k— p)-x 



1 /■ d fcdp r t (i)2) t t r (i,3) +T (M) t 

2 7 (27r) 3 L fe ' p k p k,P a k a ~P^ 1 k,p a -k a P^ 1 k,p a -k a -p 



= i(k-p)-x 



(20) 



where the and T^) are calculated from the bilinear product of the matrices S kjP and 
Tj[p respectively, with the vectors defined in (19). In taking the VEV of (20), only the term 
with the correct order survives. Then the 5 function allows us to integrate over one of the 
two (equal) momenta. Moreover, in order to handle with this expression in a cosmologi- 
cal context (as we are going to do), we assume that particles are created homogeneously 
throughout the space, so that the integral in (20) depends only on the square momentum 
k 2 of the excitation. We thus get 



r+oo 

(0| Pm (7 ? )|0) = 2vr / dkk 2 sl%, m ), 
Jo 

r+<x> 

(0\PmAv)\0) = 27r / dkk 2 Tl l f(r,, Vo ). 
Jo 



(21) 



The coefficients and T^' 4 ) are easily calculated. After a little algebra we get: 



S ( k %,Vo) = ^- {k 2 + V( V ))\A k ( v , Vo )\ 2 + \C k (v,Vo) 



2u k 



+ 



2 L 



1 

2uj k 
u k 



(k 2 + V( V ))\B k ( V , Vo )\ 2 + \D k (ri, m )\ 2 
[k 2 -V( V )) \A k ( v , Vo )\ 2 + \C k ( v , Vo )\ 2 
(k 2 -V( V ))\B k ( v , Vo )\ 2 + \D k (ri, m )\ 



(22) 



This expression can be further specified by noting that, if the space-time is flat in 7] = 7]q, 
solutions for the evolution equation (9) can be written as: 



fk(v -> Vo) = 
9k{v -> Vo) = 



/ 2^fc 



-iu k {rj-r)o) 




(23) 



with uj k = Vk 2 + m 2 (m is the mass of the excitation). Furthermore, we get a better insight 
on the expressions of energy density and pressure by writing the mode functions f k and g k 
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in terms of time-dependent Bogoliubov operators [19,23]: 

a k (r)) ( ■ P A ( \\ fofo) (■ Pa l ^ 
v / 2w fc (r/) V / V 2uj k{V) V J^, / 

Sfcfa) = - i \j^^-a k {ri)ey^{-i j dx Lv k ix^ + i\j^^p k ^) exp (^i j dx w fc (x)^ , 

(24) 

with afc and satisfying the condition \a k \ 2 — \(3 k \ 2 = 1> which follows from the Wron- 
skian condition on f k and g k . Substituting this in (22), and plugging everything into the 
expressions for the energy density and pressure (21) we get: 

r+oo 

Pm ( V ) = 2vr / dkk 2 {l + 2\p k {r ] )\ 2 )i J j k {r,), (25) 
J o 

/•+oo ^4 

p m (r?) = 2vr / dk — (l + 2\(3 k ( V )\ 2 ). (26) 

In the next section, we will use this expressions to obtain the energy density as a function 
of the geometry of the on-shell universe 



3 Energy density of the gravitationally excited string states 

In the previous section, we obtained a formal expression for the energy density and pressure 
of massive scalar modes during the evolution of the universe. At this stage we can already 
notice that, since (26) is suppressed by a factor m with respect to (25), the pressure will 
be negligible with respect to the energy density for large mass (non-relativistic) modes. 
Therefore, from now on we will assume that the pressure is exactly zero, leaving a more 
detailed treatments for future works. 

Asymptotically, the Bogoliubov coefficient f3 k (rj) can be interpreted as the number den- 
sity of particles produced via gravitational interaction from an initial vacuum state. Let us 
stress that, by setting (3 k {v) = 0> which is true for rj — > rjo, expression (25) reduces just to 
the unrenormalized zero point energy of the field. Thus we are led to consider a properly 
renormalized form for (25) as the energy density of massive excitations of the field under 
consideration, generated by gravitational interaction. We choose to regularize the energy 
density by subtracting the time-dependent zero-point energy of the oscillations, so to get 

r+oo 

Pmiv) = 47r/ dk k 2 \p k ir!)\ 2 u k iv). (27) 

J 

Such expression for the energy density has the key feature of vanishing at v = Vo (as it 
is expected since we want a vacuum solution to start with), and in general whenever the 
space-time is flat. In addition, it has the expected form of a sum over energies of each 
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mode 1 , following the interpretation of (3 as the number density of gravitationally produced 
particles, with the spectrum modified by gravitational interaction. 

We now need to evaluate fa(rj), which is in general a very difficult task. If we assume that 
the evolution of the Universe is "slow" enough, so that we are in an adiabatic regime [17], 
we can use for fa(rj) the approximate expression [20] 

= j ex P ("2« /" ^Wr/)) • (28) 

The integral can be evaluated, to the same level of approximation, with the steepest de- 
scendent method [19,20] giving: 



|A(r?)| 2 = exp 



(29) 



where _£f e ff and R e g are respectively the Hubble parameter and the scalar curvature, calcu- 
lated starting from an effective scale factor 



a cS (v) = Ja2fa)-4rV (30) 

With this expression for fa, the integral in (27) converges and can be evaluated exactly. 
After some technical manipulations we get: 



7T 



'GH^-R^irt) 2m \ 2m 



(31) 

where is the Hankel function of the first kind. For heavy massive states we can 

approximate the Hankel function with its large argument expansion [24], thus getting: 

As expected, the contribution for each massive mode is exponentially suppressed at 
low (sub-Planckian) energies. Incidentally, being at sub-Planckian energies justifies our 
quantum field theory approach to production of string modes (which, of course, have masses 
above that limit), as we stressed in the introduction. Now, as we already mentioned above, 
the multiplicity of string modes is exponentially growing [19] iV oc exp(m/T#) (T# being 
the Hagedorn temperature), and this can possibly compensate the exponential suppression 



1 Note that uJk(r/) is not exactly the energy of the single mode, since it contains also the "pump field" 
term 



7 



of each produced mode. In fact, the total energy density generated by all the stringy massive 
modes is: 

Ptotiv) = Prn(v)e*P l^J- 

m>m 3 ' ' 



(33) 



Again, this expression can be estimated by noting that a e fj(?7) ~ a(rj) for heavy enough 
states, and approximating the series as an integral, substituting m ~ Tjjx, with x > 0. We 
get: 



_3 , / r>\ 3/4 r +oo 

Ptotiv) - ^V(V-|j J d,x 5 / 2 exp 



- 1 



(34) 



Since the energy density is a scalar, turning to cosmic time is straightforward. We finally 
get: 



Ptot(t) * T* H a\t) 



H\t) 



R(t) \ V2 



T 2 



(35) 



In the next section, we will use this expression as a source in string cosmolgy equations. 



4 Cosmology with the backreacting energy density 

In the last section we have found an anaylitic expression for the energy density massive 
string states excited by the gravitational interaction with an evolving universe. Now we 
assume that this energy density backreacts on the geometry of the Universe to favor a 
smooth transition between a phase of super inflationary accelerated expansion with growing 
curvature to a (standard) phase of decelerated expansion and decreasing curvature. To see 
this, we plug the expression (35) into the pre-big bang equations. We assume that the 
extra dimension are compactified down to the Plank scale, and there is a mechanism for 
the stabilization of the modula. The string cosmology equations then read [3] 

^-6H4> + 6H 2 = ^e+ptot, (36) 

2'<p + 6H<p -<p 2 -6H - 12H 2 = 0. (37) 

(we have identified the string mass with T#) where p to t is given by (35) and we have assumed 
that matter is minimally coupled to the metric of the string frame, without any direct dilaton 
coupling, see, e.g. [25] (we will comment more on this later). Of course, these equations 
cannot be solved analytically, so we must turn to numeric. Plots for the interesting physical 
quantities are shown in Fig. 1. Numerical integration has been carried out by fixing a 
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t o'.l 0.2 0.3 0.4 0^5 0.6 0.7 

(c) (d) 

Figure 1: Plots of the dilaton energy density </> 2 (dotted line) and the squared Hubble pa- 
rameter H 2 (dashed line) (a), the Ricci curvature R (b), the energy density of the generated 
massive string states p (c), and the parametric plot of the "Shifted dilaton" energy density 
versus the Hubble parameter (d). The Hagedorn temperature has been set to 1. 

dummy initial timescale so that the maximum value of the curvature occurs at t = 0, 
and by setting suitable (small, in units of string mass) initial condition for the Hubble 
parameter and the dilaton energy density. As we can see from panel (b), the curvature 
reaches a maximum and decreases, evolving towards a phase of De Sitter evolution (i.e. 
with our notations, a constant curvature regime). It is possible to check, in addition, that 
all the approximation requirements are met. On the other hand, the dilaton energy density, 
depicted in panel (a) with the squared Hubble parameter is not under control and blows up, 
after reaching a first local maximum at the same time age of the curvature (see also plot (d)). 
This is not surprising, since we have not considered the coupling of the generated massive 
states with the dilaton and the (possible) corresponding dilaton stabilization. Motivated 
by the successful transition to a phase of decelerated expansion we found in our approach, 
we seek to extend the present analysis to the dilatonic sector in a forthcoming paper, so to 



9 



have a full-fledged graceful exit solution. 

5 Conclusions and Outlook 

Particle creation generated by gravitational backreaction seems to provide a cogent mecha- 
nism to stop superinflationary accelerated expansion and to drive the Universe to a (stan- 
dard) phase of decelerated expansion. Calculations we have presented in this paper, though 
being speculative in some of their aspects, show that, as the curvature grows, there is the 
possibility that gravitational energy excites heavy string states, thus producing some sort 
of "friction" effect. This can actually slow down and eventually stop the growth of the cur- 
vature and prevent the occurrence of the singularity. This can happen before the curvature 
reaches the string scale, ao before the universe enters in a regime in which quantum field 
theory and supergravity are no longer valid. 

However, the analysis we have performed needs to be further developed in many direc- 
tions: First of all, some of the mathematical derivations need to be posed on a more steady 
ground (for example, our derivation rely on the assumption that the metric is conformally 
flat, so that formal results valid in a flat space can be straightforwardly extended to our 
curved model). In addition, as we already noted, the particle creation mechanism has to 
include a direct coupling to the dilaton. In fact, even though the geometry is regular every- 
where in the toy model we have presented, the dilaton still blows up, indicating a breakdown 
of the validity of our approximations. Nevertheless, since the dilaton itself is a string mode, 
it is expected to couple with massive states, so high energetic dilaton modes will somehow 
excite heavy modes as well. This will probably provide a mechanism to stabilize the dilaton 
after the transition. 

Moreover, to actually enter into a phase of standard evolution, one should consider some 
sort of "reheating process" . The idea is that string states generated via gravitational back- 
reaction decays into radiation and matter, to start the standard phase. Another issue worth 
to be pursued is to have a complete dynamical system analysis. In fact, our calculations 
confirm the known result that late-time evolution of the Universe falls toward a de Sitter 
phase. It is not clear, however, if the de Sitter curvature scale is controlled by the string 
scale, or can be influenced by backreaction. In the second case, maybe one can trigger the 
mechanism (in addition to other contributions coming from a correct approach to modula 
and extra-dimension stabilization) to lower the de Sitter scale towards some realistic value. 
In this way we could possibly use the same tools that cure the pathological behaviour of 
the early Universe also to address problems of our late time Universe. 
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